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O ' Abstract 

The impact of relativistic pionic correlations and meson-exchange currents on the response 
functions for electromagnetic quasielastic electron scattering from nuclei is studied in detail. Results 
^ I in first-order perturbation theory are obtained for one-particle emission electronuclear reactions 
C ■ within the context of the relativistic Fermi gas model. Improving upon previous analyses where 
non-relativistic reductions of the currents were performed, here a fully relativistic analysis in which 
^ I both forces and currents are treated consistently is presented. Lorentz covariance is shown to play 
' a crucial role in enforcing the gauge invariance of the theory. Effects stemming uniquely from 
relativity in the pionic correlations are identified and, in particular, a comprehensive study of the 
self-energy contributions and of the currents associated with the pion is presented. First- and 
second-kind scaling for high momentum transfer is investigated. 
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1 Introduction 



For years electron scattering has provided one of the most powerful means to study the 
structure of nuclei. In particular, inclusive (e,e') processes at or near quasielastic peak 
kinematics have attracted attention in the last two decades and a number of experiments have 
been performed with the aim of disentangling the longitudinal and transverse contributions 
to the quasielastic cross section. Calculations based on a simple picture of the dynamics of 
the quasielastic peak (one-photon-exchange, impulse approximation and one-body currents) 
appeared to be successful in explaining the early (e,e') cross section data for several nuclei 
except at large energy transfer where A excitation of the nucleon contributes substantially. 
However, a simultaneous explanation of the separated longitudinal, , and transverse, , 
response functions could not be fully reached when representing the quasifree peak simply 
as an incoherent sum of elastic nucleon scattering processes. Some authors attributed this 
problem to a large quenching thought to occur in R^ |]^, ^ ^ ^, although others found 
no strong evidence for such a quenching [0, 0, || and later analysis of the world data 0] 
showed that the reduction of the longitudinal response at high momentum transfers was 
actually smaller than previously assumed. Thus at present it appears more likely that the 
"quenching" of the longitudinal response is not where the problem lies, but rather that the 
transverse response has contributions other than those mentioned above — some of these 
contributions constitute the focus of the present work. 

Different approaches to this problem have been invoked. As mentioned above, the earlier 
approaches used only one-body (nucleonic) currents, i.e., the impulse approximation (lA). 
Whereas the lA explains the electron scattering reaction mechanism around the quasielastic 
peak reasonably well, theoretical models based on one-body currents are unable to account 
for the observed strength in the dip region between the quasielastic and A peaks, where 
meson production, including via the A, and two-body currents [10| should be taken into 
account in the analysis. Several studies dealing with the issue of two-body currents in 
quasielastic electron scattering reactions have been published. Most focused on the role 
played by the meson-exchange currents (MEC), specifically the seagull and pion- in- flight 
terms. For instance, detailed results were obtained with the Fermi gas model (RFG) ||Tl|, 



12, 13] and, as well, using a continuum shell model |14, 15, 16 1. Moreover, the contribution 



of the A isobar current, i.e., the current associated with the electromagnetic excitation of 
an intermediate isobar, which decays by exchanging a virtual pion with a nucleon, has been 
also considered WA, ITB, O and shown, in most cases, to be larger than the purely nN MEC 
contributions. 

In some previous work [|l^ a relativistic analysis of the MEC has been presented (although 
there the associated correlations were neglected and thus gauge invariance was not respected); 
however, most calculations have been done basically within a non-relativistic framework. 
As a consequence, not only have non-relativistic wave functions been used, but also non- 
relativistic current operators, the latter being obtained using standard expansions which view 
both the dimensionless momentum transfer k = q/2m and the dimensionless energy transfer 
A = uj/2m as being small (m is the nucleon mass). Obviously, for high-energy conditions, the 
current operators so obtained are inadequate and new expressions are needed. This provided 
the focus for the developments in |19[ where new current operators that are exact as far as 
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the variables u and q are concerned were derived. 

As mentioned above, another difficulty of most of the non-relativistic pionic calculations 
existing in the literature relates to the breakdown of gauge invariance. In this connection, 
we recall some differences between the pionic and isobar two-body current operators. In 
a representation of the nucleus with states described only in terms of nucleons and where 
current and interaction effects are embodied in operators that contain all other hadronic 
degrees of freedom (pions, deltas, ...) these two enter somewhat differently. In lowest order 
the latter is gauge invariant by itself, whereas the former enters both in one-boson-exchange 
interactions {i.e., potentials in the non-relativistic limit) and in one-boson-exchange currents. 
Thus the two-body pionic currents must, via the continuity equation, be connected with the 
interaction effects, the so-called correlations. One might expect such pionic contributions 
to be more constrained than those arising from the isobar. But in most non-relativistic 
calculations the contribution of vr meson-exchange currents violates current conservation. 
This constitutes a fundamental flaw for the consistency of the theory and hence is one of the 
main issues addressed in this work. 

In this paper our aim is to investigate the role of pions in inclusive electron scattering from 
nuclei within the kinematical regime of the quasielastic peak. We restrict our attention to the 
study of the purely electromagnetic response functions, and i?"^, leaving for subsequent 
presentation the analysis of parity- violating responses and of the roles played in the responses 
by mesons heavier than the pion. Furthermore, here our focus is placed on the contributions 
that affect the one-particle one- hole (Ip-lh) sector of the nuclear excitations, whereas in the 
near future we shall also present results within the same framework for contributions that 
affect the 2p-2h excitations. 

This work extends the approach presented in a previous paper |T^, but introduces some 
important new elements, not only in connection with the implementation of Lorentz covari- 
ance, but also as far as the self-consistency of the theory (specifically the fulffilment of gauge 
invariance) is concerned. As in 0] we undertake this study within the context of the RFG 
model, which is well suited for maintaining the consistency in the treatment of the forces and 
currents while fullfilling the Lorentz covariance. The RFG is in fact the simplest model that 
provides a valid starting framework for describing quasielastic electron scattering and ex- 
actly respects relativity. To employ this model and to extend it using standard perturbation 
theory to account for ffist-order pionic effects amounts to considering all Feynman diagrams 
carrying one pionic line (a discussion of the relevance of higher-order pionic diagrams can be 
found in [^, ^). Accordingly, our theory is based on these diagrams and we do not make 
any non-relativistic expansion when handling the pionic current operators; herein lies the 
main difference with previous work ||T2| , p!9| , where some kind of expansion in powers of 
1/m was introduced. In the present study we present a fully relativistic analysis including 
the pionic correlation effects, embodied in the so-called self-energy and exchange diagrams, 
as well as the MEG contributions (seagull and pion-in-fiight terms). The analysis of the A 
current, within the same framework, will be presented in a forthcoming publication. 

With regard to gauge invariance, in contrast to previous (non-relativistic) analyses, our 
findings strongly support the need for a fully relativistic treatment in order to preserve this 
symmetry. In addition, relativity is clearly manifested in new effects which appear at the 
level of the currents and/or propagators and which are absent in the non-relativistic approxi- 
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mations. The relativistic analysis of the self-energy (s.e.) term that presents some important 
differences compared with the non- relativistic approach is particularly significant. This issue 
will be treated in detail in forthcoming work |23| and so here we only mention the fact that 
two different contributions in the s.e. pionic correlation term can be identified. One stems 
from the renormalization of the energy and the other has its origin in the lower components 
of the Dirac wave functions. Whereas the first term can be shown to correspond to the one 
already embodied in the s.e. diagram within the non-relativistic approach |T^, the second 
one has no counterpart in the non-relativistic analysis. However, both contributions are nec- 
essary, not only to set up a consistent and well-defined self-energy particle-hole (ph) current 
matrix element, but also to fulfill the continuity equation expressing the gauge invariance of 
the theory. 

This paper is organized as follows: in Section ^ we present the general formalism and 
the basic ingredients that enter in analyses of inclusive quasielastic electron scattering. We 
define the longitudinal and transverse electromagnetic inclusive response functions in terms 
of the hadronic tensor W'^'^. This in turn can be constructed through the ph current matrix 
elements or, equivalently, as the imaginary part of the polarization propagator. Section |^ 
deals with the one-pion-exchange contributions to the nuclear responses. We start from the 
Feynman diagrams corresponding to the free-space MEG and pion induced nucleon-nucleon 
correlation currents. We prove, at this level, the gauge invariance of the theory. From these 
currents we next set up the ph matrix elements that enter in our RFG-based computation. 
The case of the s.e. diagram is discussed in detail as the corresponding ph current matrix 
element is divergent. We treat this term using the polarization propagator approach, thus 
obtaining a finite response. In p3[ we shall prove, however, that it is possible to derive 
a new "renormalized" expression for the self-energy current matrix element which is finite, 
consistent with the analysis performed with the polarization propagator technique and, when 
taken together with the remaining pionic currents and the one-body current, satisfies current 
conservation. In Section^ we discuss the effects introduced by the MEG and pion correlations 
at the level of the two response functions, and R^. We study the magnitude of these 
effects as a function of the transfer momentum q, and, as well, their dependence upon the 
Fermi momentum kp. Such dependences are clearly relevant in connection with analyses 
of scaling in nuclei |]2^, both of first and second kinds, and, as we shall see, scaling of 
the former kind is attained at high q, whereas scaling of the latter kind is clearly broken 
through dependence roughly on the Fermi momentum squared for all values of q. The issue 
of the impact on the responses of the dynamical pion propagator, versus the standard static 
propagator, is also discussed. Finally, in Sections ^ and ^ we summarize our results and 
draw our conclusions. 



2 General Formalism 



In this section we briefly summarize the general formalism involved in the description of 
(e,e') processes for quasielastic kinematics. We only discuss those aspects that are of special 
relevance to the analysis that follows and the reader who is interested in more details will 
benefit from reading |2|, |6|, 0, |28|, |9 



We limit our attention to the Plane Wave Born 
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Approximation (PWBA) in which the electron is described as a plane wave and interacts with 
the nuclear target via the exchange of a virtual photon. The variables involved in the process 
are Kj^ = [e, k) and K'f- = {e', k'), the initial and scattered electron four-momenta, and Pf^ = 
(Mj, 0) and Pj = {Ef,pf), the initial and final hadronic four-momenta, respectively. The 
four-momentum transferred by the virtual photon is Q'^ = {K^. — K'^^ = [Pf — PiY = (u^, q)- 
Assuming Lorentz invariance and parity conservation, the inclusive (e,e') cross section then 
reads 



da 1o? 



T. 



where a is the fine structure constant, aM the Mott cross section and Q'^ the scattered electron 
solid angle. The terms rjfj,iy and W'^'^ represent the leptonic and hadronic tensors, respec- 
tively. The kinematic factors vl and vt are evaluated from the leptonic tensor with standard 



techniques (see p5[ for explicit expressions), whereas the longitudinal and transverse (with 
respect to the momentum transfer q) response functions and R^ are constructed directly 
as components of the hadronic tensor W^'^ according to 



q q^ 



R:^(g,cj) = W" + W'''. (3) 

Note that if gauge invariance is fulfilled, implying that W'^^ = W^'^ = {uj/q)W'^^ and 
W^^ = {u!/q)'^W^^, then R^ is simply the time component of the hadronic tensor, namely 
W^^. Hence R^ is determined by the charge distribution, whereas R^ reflects the current 
distribution of the nuclear target. 

The hadronic tensor and consequently the response functions derived from it embody 
the entire dependence on the nuclear structure, specifically on the charge and current dis- 
tributions in nuclei, and accordingly these provide the prime focus in analyses of electron 
scattering. There are various options on how to proceed in performing such analyses (see, for 
example, depending on the specific problem under consideration and on the approx- 

imations to be made. In what follows we recall two common expressions for the hadronic 
tensor W^" and comment briefly on their applications. 

First, the hadronic tensor can be defined according to 

= E E(/l'^1^)*(/l^1^)^(^^ + uj-Ef), (4) 

i f 

where J'^ represents the nuclear many-body current operator, the nuclear states \i) and |/) 
are exact eigenstates of the nuclear Hamiltonian with definite four-momentum, and the sum 
with a bar means average over initial states. This form is very general and includes all 
possible final states that can be reached through the action of the current operator on 
the exact ground state. In our perturbative approach we shall use eigenstates of the free 
Hamiltonian Hq (which describes the free RFG) and include correlations among nucleons in 
the current mediated by the exchange of pions (MEG). This current of course allows one to 
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reach both the p-h and the 2p-2h sectors in the Hilbert space of Hq. In the present work, 
however, we shall restrict our attention to the former. 

A different option for evaluating the nuclear responses exploits the polarization propa- 
gator n^'^ (also referred to as the current-current correlation function). The latter can be 
expressed in terms of the full propagator, G, of the nuclear many-body system, since closure 
can be used to carry out the sum over the final states in eq. (^. Then one has for the 
hadronic tensor 

W^"" = --Im (W) = --Im ( y + EArH)] . (5) 

^ ^ \i J 

A possible advantage of this approach relates to the existence of a well-defined set of rules 
(the relativistic Feynman diagrams) which allows one to compute U'^'^ perturbatively . 



Obviously both procedures are equivalent and hence the observables calculated using the 
expressions for the hadronic tensor given by eqs. (^) or (|^) should be the same. However, 
notice that eq. (^ is less suitable for dealing with situations where the nuclear current matrix 
element {f\J^\i) is divergent. In this case it is advisable to compute the responses via the 
polarization propagator. 

Finally, we remark that gauge invariance must be fulfilled both at the level of the nuclear 
current matrix elements and at the level of the hadronic tensor and/or the polarization prop- 
agator. A consequence is that the electromagnetic continuity equation should be satisfied. 
In other words in momentum space all of the expressions Q^iW^" and Q^jy^" 

should vanish. 



3 Pion exchange contributions 

In this Section we shall discuss the hadronic tensor and response functions in the framework of 
the RFG model, i.e. for nucleons moving freely inside the nucleus with relativistic kinematics, 
accounting for the effects introduced by pions in first-order perturbation theory (one-pion- 
exchange). We shall employ pseudovector (PV) coupling for the pion-nucleon interactions. 



3.1 Feynman diagrams and two-body currents 

We start by displaying in Fig. 1 the Feynman diagrams corresponding to the free-space 
two-body currents with only one-pion-exchange. Using the Feynman rules one writes down 
the corresponding expressions for the two-particle current matrix elements. Following stan- 
dard terminology, diagrams (a), (b) and (c) correspond to the usual meson-exchange currents 
where the electromagnetic field is coupled with the pionic current: in particular, diagrams 
(a) and (b) are referred to as the seagull (contact) term and (c) as the pion-in-flight term, 
respectively. Diagrams (d)-(g) contain an intermediate virtual nucleon and give rise to the 
so-called correlation currents. In some work it has been argued that the contribution of 



these "correlation" diagrams is already contained in the nuclear wave functions as solutions 
of the Schrodinger equation. However, within the RFG plus first-order perturbation theory 
framework the unperturbed basis does not have any such contributions and so within this 
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Figure 1: Free-space two-body current diagrams: Seagull (a,b), pion- in-flight (c) and corre- 
lation (d-g) contributions. 

approach they must be included explicitly — in fact, as we shall see, their inclusion is crucial 
for fulfilling gauge invariance. 

The relativistic expressions for the free-space two-particle currents are (isospin summa- 
tions are understood) 

• Seagull or contact: 



fm 



\VmJ 
+ (1^2) 

Pion-in-flight: 



F 

i^3abU{p'l)Ta'y5 ^iu{pi)—^ ^M(P2)^fe757M«(P2) 

-"■1 ~ ^TT 



(6) 



fm 



X M(pi)Ta75 /^iM(Pi)«(P2)t675 ^2«(P2) 



2)11 



(7) 



Correlation: 
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+ (1^2). (8) 

In the above, V is the volume enclosing the system, is the mass of the pion, / the 
pion-nucleon coupling constant and Ep = + the on-shell energy of a nucleon with 
momentum p; the four-momenta — indicated by capital letters — are defined in Fig. 1 and 
F^ and are the electromagnetic isovector nucleon and pion form factors, respectively. 
Furthermore, Sf{P) is the nucleon propagator and T'^{Q) the electromagnetic nucleon vertex, 
i.e., 

Sf{P) = (9) 

r'^(Q) = Fi7^ + -i-Fsa^-Q, , (10) 

Fi and F2 being the Dirac and Pauli form factors: for these we use the Galster parame- 
terization IQ. Finally, the spinors (for brevity we denote u(p, Sp) by m(p)) are normalized 
according to the Bjorken and Drell convention |^4[ and isospinors are not explicitly indicated. 
The seagull and pion-in-flight currents shown above coincide with the expressions given 



by Van Orden and Donnelly |11| if account is taken for the different conventions used for 



the gamma matrix 75 and for the metric. Concerning the correlation current, note that. 



at variance with |]TT|, it embodies both the positive and negative energy components of the 
nucleon propagator. 

An important point to be stressed is that the relativistic seagull, pion-in-flight and cor- 
relation currents satisfy current conservation, i.e. Q^J^ = 0, provided some assumptions 
are made for the form factors involved in the various currents. To prove this statement we 
start by evaluating the contraction of the four- momentum transfer Q'^ with the correlation 
current . It can be written as 



with JUa given by 

Ma = m(p^) [r,75 ^iSf{P2 + Q) ^Fi + F, QSf{P^ - Q)Tai, ^1] ^(ps) , (12) 

where we have used the relation Q'^T^^^Q) = Fi{Q) Q. After some algebra, involving the 
nucleon propagator and the Dirac spinors, M.a can be further simplified leading to the form 

Ma = U(P2) Ms J^lFl - F^Tal^ ^i] n(p2) = ui^'^) [fa, F^] 75 ^ Mp^) • (13) 

In order to evaluate the commutator [r^, F]\ we decompose the nucleon form factor into its 
isoscalar and isovector pieces, Fi = ^ [f^ + F^r^^ ■ Then 

[Ta,F^] = -iF^e3abn : (14) 
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which entails the automatic conservation of the tc^ exchange current (a=3). Using eq. (|T^) 
we can recast as follows 

Ma = -iF^e3abU{P2)ni5 ^MP2) ] (15) 

hence the divergence of the two-body correlation current matrix element can finally be 
written as 

Q^^i^ = 2m — — , _ _ ^ ie3abu{p'i)ral5uiPi^ ^ 



VmJ ^E^,E^,E^,E^,^ ' Kl - ml 

X u{p'^)n^^{Q + 2m)u{p2) + (1^2). (16) 

The divergence of the seagull and pion-in-fiight two-body current matrix elements can 
also be calculated in a straightforward way. The final result reads 

X «(P2)rj75 ^«(P2) + (1 " 2) (17) 



Vm^) ^Ep,Ep,Ep,^Ep,^ " {Kf - ml){Kl - m^ 

X u{p[)Ta'j5u{pi)uiP2)n'j5u{p2) ■ (18) 

Then, by summing up the contributions given by the correlation (eq. (|16])) and seagull 
(eq. (p!?!) ) currents and writing the four-momentum transfer as = {Ki + K2)^, we finally 
obtain 

X M(p'i)ra75M(pi)M(p2)rb75M(p2) , (19) 

which cancels exactly the contribution of pion-in-fiight current in eq. ([T8|) provided the 
electromagnetic pion form factor is chosen to be F.,^ = , which we shall assume in this 
paper. 

We have thus proven that the sum of the relativistic two-body matrix elements of the 
MEC and of the correlation currents satisfies the continuity equation if the same electro- 
magnetic form factors enter in all of the currents: namely one has Q^{j^ + + j^) = 0. 

3.2 Particle-hole matrix elements 

In this paper we focus on the one-particle emission reactions within the framework of the 
RFC Therefore, one has to evaluate the matrix element of the two-body current operators 
taken between the Fermi sphere and a particle-hole state, namely 

{ph''\j,\F) = E E [{Pk\j,\hk) - {pk\j,\kh)] , (20) 

Sk,tk l<.<kF 
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Figure 2: Feynman diagrams representing the seagull (a and b), pion-in-flight (c), vertex (d 
and e) and self-energy of the hole (f) and particle (g) particle- hole matrix elements. 



where the sum X]k<A;F the thermodynamic limit becomes an integral over the momentum 
in the range < k < kp, with kp the Fermi momentum, and over the angular variables 
9k,4>k- The first and second terms in eq. (^) represent the direct and exchange contribution 
to the matrix element, respectively. It is well-known that the direct term vanishes for the 
MEC (seagull and pion- in- flight) currents (see [0) because of both spin-isospin traces (in 
spin-isospin saturated systems) and the presence of a pion carrying zero momentum. It is 
easily verified that, for the same reasons, the direct matrix element of the correlation current 
also does not contribute to the responses. 

Hence only the exchange term remains in the particle-hole matrix elements for the two- 
body current operators. The corresponding Feynman diagrams are displayed in Fig. 2. 
Diagrams (a-b) and (c) correspond to the seagull (contact) and pion-in-fiight terms, respec- 
tively. Diagrams (d-e) represent instead the correlation contributions. Here we distinguish 
the self-energy insertions on the nucleonic lines (f,g) and the exchange of a pion between a 
particle and a hole line (d,e). 

Starting from the general expressions in eqs. (|,0,||) we can now evaluate the particle-hole 
matrix elements. Carrying out explicitly the sums over the spin, s^, and isospin, tk, variables 
and after some laborious, but straightforward, algebra we finally get 



Seagull 
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V^ml^EpEh 

X E -w(p)r„.. I + _ I .(h) (21) 



Pion-in-flight 



m {Q + 2H - 2K)^ 



ab 



Correlation 



X [r,75(/'- ^)5^(// + Q)r^(Q) + T,{Q)Sf{K - Q)r,75(f- ^)] 
+ [r,75(^- ^)S^(ir + Q)r^(g)+r^(g)^^(P-Q)r„75(^- /f)] 

X (^ + ^K75^^^-^}«(h) . (23) 

The first and last terms in eq. ( p3D represent the so-called self-energy (s.e.) matrix elements 
(diagrams (g) and (f) in Fig. 2), whereas the second and third terms correspond to the vertex 
correlations (v.c, sometimes referred to as "exchange") matrix elements (diagrams (e) and 
(d)). We examine them separately, writing accordingly 

{ph-'\f,\F) = {ph''\f,\F)s.^. + {ph-'\f^\F)^.^. . (24) 

The self-energy ph matrix element can be split into two parts, 7^^ and Til^. The former 
corresponds to the diagram with the pion inserted in the particle line (Fig. 2g), whereas the 
latter describes the diagram with the pion inserted in the hole line (Fig. 2f). They are 



K= ^^^ ^^jl ^ Y.^u{v)^'iQ)SF{H) P -m)u{\.) , (26) 



where the isospin trace has been performed. 

Similarly, the vertex ph matrix element splits into two terms and B'^ representing the 
forward- and backward-going contributions, respectively (Fig. 2d and 2e). They are 
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(27) 

(28) 

Finally, splitting also the electromagnetic nucleon operator into its isoscalar and 
isovector parts, one obtains the isoscalar and isovector contributions to the self-energy and 
vertex ph matrix elements. The final results can be cast in the form 



(29) 

and 



(30) 

for the self-energy matrix elements and 



3m/2 



X E £^(P)75(^- mSF{K + g)F^(")(g)75 ^. uih) (31) 

5 ,,, ^2 



^ ■3"'/ 



V'ml^E^E^ 



X E |^^(P) rP ^K\T ^2 ^^"^'^'^^)^^^^ - ^)75(/'- ^)«(h) , (32) 



for the isoscalar and 
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k - Hy - ml 



(33) 



(34) 

for the isovector vertex ph matrix elements. Interestingly, as far as isospin traces are con- 
cerned, in the vertex matrix element the isoscalar /isovector ratio is -3, whereas in the self- 
energy case ■t .s 1 a Note that the MEC (pion-m-flight aud seagull) ph matrix elements 
are purely isovector, whereas the vertex and self-energy correlations get both isoscalar and 
isovector contributions. 

Note that the v.c. and s.e. ph matrix elements involve the nucleon propagator Sf{P) 
which in some situations may imply the occurrence of singularities. In the case of the 
vertex diagrams, the four-momenta appearing in the propagators are K + Q and K — Q 
for the forward- (Fig. 2d) and backward-going (Fig. 2e) contributions, respectively, and 
an integration over k should be done. For q > 2kF (no Pauli blocking) it can be proven 
(see appendix B) that only the forward diagram contains a pole, i.e., a value of the inner 
momentum k exists such that the nucleon carrying a four-momentum -|- Q is on-shell. 
In this situation the forward vertex ph matrix element is evaluated by taking the principal 
value in the integral over cosOk. In the case of the backward-going diagram the nucleon 
propagator Sf{K — Q) has no singularity for the kinematics in which we are interested. 

The case of the self-energy diagrams is different. Here the ph matrix element diverges, 
as the nucleons described by the Feynman propagators Sp{P) and Sp{H) are forced to be 
on-shell by energy and momentum conservation. Accordingly this term should in principle 



be renormalized. In |£3| we present a detailed analysis of the different aspects of the impact 
of the self-energy term on the nuclear responses within the relativistic framework. A compar- 
ison with the non-relativistic reduction is also carried out there. Here we confine ourselves to 
computing the relativistic response functions through the hadronic tensor W^". As stated in 
Section 0, W'^'^ can be evaluated using the current matrix elements or the polarization prop- 
agator. In the first instance, however, the presence of divergent ph matrix elements may lead 
to wrong responses or inconsistencies in the theory. In this paper, to avoid this happening, 
we calculate the s.e. contribution to the response with the polarization propagator. 

An important point here is the fact that the hadronic tensor thus obtained, partly through 
the ph matrix elements and partly through the polarization propagator, is gauge invariant. 
This may be somewhat surprising because we show in appendix A that current conservation 
is already obtained at the level of the MEC and correlation ph matrix elements: hence the 
one-body current ph matrix element also has to be conserved. However, this occurs only in 



^The latter result stems from the relation T3 + ie^abTaTt — — T3; however we prefer to leave the isospin 
structure of the isovector exchange as in (js^,^) since it makes more transaparent the self-energy-exchange 
cancellation in the continuity equation, as shown in appendix A. 
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zeroth order of perturbation theory. To be dealt with properly, clearly the situation requires 
the renormalization of the ph energies and Dirac spinors. Only then does it become possible 
to set up a renormalized s.e. current which leads to a hadronic tensor that coincides with 



the one obtained here through the polarization propagator [23 



4 Nuclear hadronic tensor and response functions 

In this Section we analyze the response functions for one-particle emission reactions within 
the RFG model. As discussed in the previous Section, the ph matrix elements corresponding 
to the different pionic diagrams are all well-defined except for the self-energy term which 
diverges. Therefore, on the one hand, in the case of the one-body, MEG and vertex correlation 
diagrams we evaluate the hadronic tensor starting from the current ph matrix elements 
(Fig. 2). On the other hand, for the self-energy diagrams we calculate the polarization 
propagator U'^'^. 

The analysis of the nuclear hadronic tensor set up with the ph matrix elements has been 
presented in detail, within the RFG model, in previous papers |]T2|, |20|. Hence, here we 



simply summarize the results needed for later discussions. The hadronic tensor that arises 
from the interference of the single- nucleon current, with the one-pion-exchange current 
jj^, with a = S (seagull), P (pion-in-flight) and v.c. (vertex), is for the RFG model with 
Z = N 



= ^T^ / ^^^(^ + ^h) / d<p, Yl 2Re [{ph-'\jUFr{ph''\j:\F)] , (35) 
QTTKpq Jho JO ^ 



where {ph \F) = .J" u{p)r^u{h) is the single-nucleon ph matrix element with 

V P h 



F'^ the electromagnetic nucleon current from eq. (|TOD and {ph ^\j^\F) is the ph matrix 



element for the seagull, pion- in- flight or vertex current as given in eqs. (^1]), (12^) and 
p^), respectively. 

Note that in eq. (|35D the integral over hole polar angle, cos 6h, has been performed 
explicitly by exploiting the energy conserving 5-function. This fixes the minimum momentum 
of the hole according to 




hQ = mJeQ — l, Eq = Max iep — 2\, K\ 1 -\ A> , (36) 



where the usual dimensionless variables 



^ ^ \Q^\ Q Ep h 

^ = 7^> ^ = r^' ^ = ^^ = — ' v = — (37) 

2m 4m^ 2m m m, 



have been introduced and Ep = \Jk\ + is the Fermi energy. Moreover, the hole three- 
momentum 

h = (sin6'o COS0/1, sin 6*0 sin cos6'o) , (38) 
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(b) (c) 



Figure 3: Feynman diagrams of the free (a) and first-order self-energy (b and c) polarization 
propagator. 

involved in the hadronic tensor, must be evaluated for the following specific value of the 
polar angle 

Xe — T 

cos ^0 = • (39) 

The hadronic tensor, as was the case for the current, can be also split into isoscalar and 
isovector parts, since there is no interference between the two isospin channels. 

An important issue relates to the form factor of the ttNN vertex, F^, which incorporates 
some aspects of the short-range physics affecting the pionic correlation. In all of the above 
expressions F^ has not been explicitly indicated for sake of simplicity. Here we recall that its 
inclusion in the ph current matrix elements is not without consequences in connection with 
gauge invariance. In fact, in this case, the results quoted in appendix A are no longer valid 



unless new terms are added to the MFC (see [10, 20, 0). Lacking a fundamental theory for 



Fjr, in the calculations reported in the next Section we use the phenomenological expression 

/\2 _ ™2 

with A = 1.3 GeV. As long as the dependence upon A is not too strong the gauge invariance 
of the theory should not be too badly affected. Within a non-relativistic approach for the 
pion currents, a detailed discussion on the breakdown of the gauge invariance due to eq. (M), 



and on the dependence of the responses upon the cutoff value can be found in pO|, pi| . 

Now, as already discussed, a crucial point to be emphasized is that the self-energy ph 
matrix element is divergent. Hence it should be renormalized to be used in the evaluation of 



the hadronic tensor. Referring the reader to P3[ for a thorough discussion on this subject 



here we show how the hadronic tensor is built by evaluating the polarization propagator 
corresponding to the two self-energy diagrams shown in Fig. 3. 
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From the Feynman diagrams with the self-energy E inserted in the particle or in the hole 
lines (Fig. 3), the polarization propagator is 



dhod^h 



Tr {T''iQ)S{H)J:iH)SiH)T''i-Q)SiH + Q) 



+ r^'{Q)S{H)r''{-Q)S{H + Q)J:{H + Q)S{H + Q)} . (41) 
In the above the nucleon propagator within the RFG framework reads 



S{K) 



K — m + ie 



9{k - kp) 



+ 



d{kF - k) 



l<[2 _ ^2 _|_ _ ^2 _ j^^U 



(42) 



Concerning the self-energy one has 



S(P) 



(P - K)^ - ml 



(43) 



In the following we illustrate the calculation by treating in detail the case of the self- 
energy in the hole line. The analysis of the diagram with the s.e. insertion on the particle 
line is similar and then we shall just quote the final result. In appendix C we discuss the 
more general case with any number of self-energy insertions in the particle and/or in the 
hole lines. 

Using the nucleon propagator in eq. (^2]), we can write the polarization propagator with 
the s.e. on the hole line as follows 



n'^'^(Q) 



dhnd^h 



(2vr)4 
9{h - kp) 



eikp - h) 



+ 



[H"^ — m? + ieY {H"^ — m? — iehoY 

h + q|-A;ir) ^ 6'(fci7 -|h + q 



+ 



Then, exploiting the identity IT! 



(44) 



1 



{H"^ — m? — ie) 



d 

da 



1 



0=0 



a — — le 



(45) 



and subtracting the (infinite) vacuum polarization (see appendix C) we get 



n^'^(g) 



d 

da 



Tr 



dhnd^h 



X 



d{kF - h) 



(2vr) 



T^{Q){P + m)E(i/)(^ + m)T''{-Q){:p+ ^ + m) 



H"^ — a — m? — ieho 



9{\h + q\-kF) 
{H + Qf -m? + ie 



(46) 
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Now the integral over the hole energy ho can be computed. For this purpose we rewrite 
the denominators of the hole and particle propagators in the form 

-a-m^ - ieho = - m'^ - ieho = {ho + - ie) {ho -E^- ie) (47) 
{H + Qf-m^ + ie = {ho + qo + Ep - ie){ho + qo - Ep + ie) , (48) 

where we have introduced a mass 



m' = m'{a) = V a + 



(49) 



and an energy 



m 
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(50) 



that depend upon the parameter a and we have used the relations p = h+q, Ep = y/p'^ + m'^. 
Therefore the poles are located at the points 



ho = ^E^^ + ie 

ho = ±{Ep - ie) - go 



(51) 
(52) 



and the integral over ho can be computed by closing the integration path in the lower half- 
plane, so that only the pole at ho = Ep—ie—qo contributes. Then the polarization propagator 
is found to be 



n^"(g) 
d 

da 



/d^h 
—-T^{Q){P + m)T.{H){P + m)r^(-g) {p+Q + m) 



a=0 



X ■ 



(27r)3 

1 e{kF-h)d{\h + c^\-kF) 



ho=Ep-qo 



2Ep {Ep - qo)^ - h2 - m'2 - te{Ep - go) 
The hadronic tensor then follows from formula (|5|) and reads 



(53) 



= --im [n'^^(g)] 

vr 

Stt^Z d 



kp da 



-—T>^{Q){P + m)J:{H){P + m)r (-g) {p+Q + m) 



X ■ 



2E, 



-9{kp-h)9{\h + q\-kF) 



HEp - go) 



5{Ep-qo-Ei) , 



'ho=Ep-qo 



(54) 



where we have exploited the inequality Ep—qo > and the relation 6{x'^—xl) = 6{x—xo)/2xo 
for xq > 0. 
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Finally, using the energy-conserving delta-function and setting pq = Ep, we obtain 



kp da 



/d^h 
—-V^{Q){P + m)m)iP + m)V^i-Q) (f + ^)l 

a=U V^^'' 



2E, 



{2n 

-e{kF-h)e{\h + q\-kF) 



ho=E^{a) 

5{Ep-q,-E'^{a)) , (55) 



an expression that coincides with the one given by eq. (|93|) in appendix C. 

Using a similar analysis, the hadronic tensor with the self-energy inserted in the particle 
line turns out to read (see appendix C) 



Stt^Z d 



k% dp 



(3=0 
X 



r d^h 

Tr j ^r'^(g)(^ + m)r(-g)(/' + m)S(P) (f + m) 
1^^^0{kF - h)e{\h + q| - kF)^5{E'^{f5) - go - E^) 



(56) 



with -E'p(/5) = a/p^ -|- m'^ and m' = m'(/3) = ^/^'^Trr? ■ 

In the above expressions, after the derivatives with respect to the parameters a and /3 
are taken, the integral over the hole polar angle cos 6h can be performed analytically by 
exploiting the 5-function. Hence the s.e. contribution to the hadronic tensor can finally 
be expressed as a double integral. Since the self-energy E involves a triple integral, the 
contribution to hadronic tensor turns out to be a 5-dimensional integral, to be carried out 
numerically. 



5 Results 

In this Section we report the numerical results obtained for the pionic MEC (pion-in-fiight 
and seagull) and for the correlation (vertex and self-energy) contributions to the quasielastic 
peak in the Ip-lh sector. The calculation is fully relativistic. We have taken Z = N = 20 
and set kp = 237 MeV/c, which is representative of nuclei in the vicinity of ^°Ca. 



The 5-dimensional integrations of the MEC and correlation responses implicit in eq. (|35| ) 
have been performed numerically. The reliability of the numerical procedure has been proven 
by checking that the free RFG responses coincide with their analytical expressions (see. 



e.g., [26]). Moreover the transverse MEC responses have been compared with the non- 



relativistic calculation developed in where the seagull ph matrix element is evaluated 



analytically, while the pion-in-fiight contribution is reduced to a one- dimensional integral: 
both calculations give the same results for q and kp small. 

5.1 MEC 

With regard to the MEC we have found the following: 

• They are almost irrelevant in the longitudinal channel, whereas they are small but not 
insignificant in the transverse one. This outcome emerges from inspection of Figs. 4 and 



17 



50 100 150 200 250 200 300 400 500 600 

a;(MeV) a;(MeV) 




1000 1200 1400 1600 1900 2100 2300 2500 
w(MeV) a;(MeV) 



Figure 4: Longitudinal response versus uj. Solid: free; dashed: MEC contribution. Here and 
in all of the figures to follow the nucleus is '"'Ca with fci?=237 MeV/c and the units are 10~^ 
MeV-i at g=0.5 GeV/c (panel a), IQ-^ MeV"^ at g=l GeV/c (panel b), 10"=^ MeV-^ at 
g=2 GeV/c (panel c) and 10"^ MeV"^ at GeV/c (panel d). 




Figure 5: Transverse response versus u) at g=0.5 (a), 1 (b), 2 (c) and 3 (d) GeV/c. Solid: 
free; dashed: MEC contribution. 



18 



5, where the MEC contribution (dashed curves) to the L and T responses is displayed 
for various vahics of q together with the free responses (sohd). While in the L channel 
the MEC are hardly visible, in the T channel they contribute somewhat more, although 
the contributions there typically only amount to about 5-10%, depending upon q and 
uj (see later) . This dominance of the spin interaction between the photon and the MEC 
prompts us to concentrate on the transverse channel. 
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Figure 6: MEC contribution to transverse response versus uj at ^=0.5 (a), 1 (b), 2 (c) and 
3 (d) GeV/c. Sohd: pion- in- flight; dashed: seagull; dot-dashed: total MEC. 




Figure 7: The ratio MEC/free in the transverse channel plotted versus ^0 at kp— '2.2)7 MeV/c 
for various values of q (in GeV/c) in panel a and for various values oikp (in MeV/c) at q—1 
GeV/c in panel b. 
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In Fig. 6 we display the separate pion-in-flight and seagull contribution to for 
various values of q. It appears that the seagull term (dashed lines) is always larger 
than the pion-in-flight term (solid lines), this dominance increasing with q and again 
reflecting the spin nature of the photon-MEC interaction. Moreover whereas the pion- 
in-fiight term is always negative, the seagull changes sign with u;, inducing a (mild) 
softening of the response. 

In Fig. 7a we study the evolution with q of the MEC: it clearly appears that their 
relative contribution to i?^ decreases with g, but does not vanish for large values 
of q. This point is best illustrated by displaying the response versus the scaling 



variable 24 



* - ^ , (57) 

Sp — i 



using eq. (0). Indeed the range — 1 < ^ 1 "will then be common to all responses, 
no matter what is the value of q. In the figure we plot the ratio = R^ec/ ^Jree 
for q=0.5, 1, 1.5, 2, 2.5 and 3 GeV/c: we see the relative MEC contribution decrease 
in going from 0.5 to 1 GeV/c, but then it rapidly saturates at or slightly above q=l 
GeV/c, where its value stabilizes, typically around 10%. Thus, one can see from these 
results that at momentum transfers above 1 GeV/c scaling of the first kind is 
satisfied for the MEC contributions considered in this work. Note also that, for high 
q the MEC almost vanish in the vicinity of the QEP (?/'=0). 

In Fig. 7b we investigate the kp dependence of the MEC contribution. Importantly, 
the latter is seen to grow with kp, in contrast with the free response which decreases 
as kp^. This point is again best studied by displaying the response as a function of ip, 
since the response region broadens with kp. In Fig. 7b the ratio is plotted versus 
ip at q=l GeV/c for kp varying from 200 to 300 MeV/c. It clearly appears that the 
relative MEC contribution grows with kp (attaining a value of about 20% at kp =300 
MeV/c). 



Indeed, our earlier studies of the current matrix elements |19l can be used in analyzing 
the functional dependence on kp. There a "semi-relativistic" approach was followed 
in which, once the leading power dependence on the dimensionless variable h/m was 
identified, where h = |h| is the three-momentum of any particle below the Fermi 
surface, all higher-order terms in h/m are neglected, since h/m is bounded by ?7i;' = 
kp/m and the latter is generally quite small. For the matrix elements this was shown 
in |T9| to be an excellent approximation. Here we see a detailed verification of those 
ideas. In fact, a full analysis at very small Fermi momentum shows that the seagull 
contribution has a leading power of kp, while the pion- in- flight contribution begins 
as kp — the former because it involves one- and two-body spin currents, the latter 
because it involves convection currents and therefore one higher power of momentum 
in both the one- and two-body contributions. Thus, if we divide the fully-relativistic 
results obtained in this work by these respective factors of kp with n = 2, 4 for seagull 
and pion-in-fiight, respectively, we should see very httle residual dependence on kp, at 
least at very low Fermi momentum. Indeed, even for values of the Fermi momentum 
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that are typical of nuclei across the periodic table this is still roughly the case, as can 
be seen by examining Fig. 8, panels a and b. The relatively small remaining spread 
with varying kp reflects the quality of the semi-relativistic treatment in which only 
the leading kp dependence is taken into account. Of course, in the present work no 
approximation is made at all and it is not necessary to use the rough power dependence 
except as an orientation to assess the degree to which the two-body MEC effects depend 
on the density in a manner that is different from the behaviour seen in the one-body 
currents. 

It is also interesting to note at this juncture that, if we attempt to find the best choice 
for the power n in the two cases, with n free to be non-integer, specifically for nuclei in 
the range 200 < kp < 260 (i.e., excluding only the lightest nuclei), then the respective 
powers turn out to be somewhat lower than 2 and 4, namely, about 1.5 and 2.5. 

Given that once the leading powers of kp are removed the integrals remaining in 
the two cases are of roughly similar magnitude, one might expect the pion-in-flight 
contribution to be much smaller than the seagull contribution, at least at high q where 
such arguments can be made relatively easily (at low q the pion propagators, etc. 
make the argument less secure). Indeed, having powers of kp and kp, as discussed 
above, would imply that at high q typically the seagull will win by roughly a factor of 
{m/kpY — 16, and this is borne out in Fig. 6. 

Thus, overall the seagull contribution is dominant and so an overall leading-fc^ depen- 
dence of about k'^p is appropriate. Since the RFG has a factor of k~^^ this means that 
the relative dependence — these two-body Ip-lh pionic contributions compared with 
the RFG — goes approximately as k^p. That is, when studying the scaling behavior 
of the second kind, these two-body MEC processes violate the second-kind scal- 
ing by roughly three powers of kp. Clearly the effect is a rapid function of the 
Fermi momentum (or equivalently, of the density): for example, if one considers the 
cases ^H/^He/heavy nuclei with Fermi momenta of approximately 55/200/260 MeV/c, 
respectively, then the Ip-lh MEC contributions amount to 0.1/5/10% of the total 
transverse response, respectively (normalizing to 10% for the heavy nucleus case — see 
the first bullet). 

In Fig. 9 we compare the MEC contribution to the transverse response with the non- 



relativistic calculation of [|T5l for q = 500 MeV/c. For this comparison we use T^^ = 1 
and the static pion propagator in the relativistic calculation. The effect of static 
versus dynamic pion propagator will be discussed in the next item. From Fig. 9 we 
see that, apart from the difference stemming from the relativistic kinematics which 
shrinks the response domain, the relativistic responses are about 30% smaller than the 
non-relativistic ones, indicating that relativity plays an important role even for not so 
high g- values. 



Finally the impact on the responses of the relativistic propagator A,r(-^) = {K 



2 



2^ 



as compared to the static one A^"'^'^(k) = — (k^ + m^)~^, which is commonly 
used in non-relativistic calculations, is explored. In Fig. 10 the pion-in-flight, seagull 
and total MEC contributions to are evaluated for g=0.5 and 2 GeV/c using both 
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Figure 8: The ratio R'^/kp of the seagull (a), pion-in-flight (b), vertex correlations (c) and 
self-energy (d) contributions to the transverse response and kp (with a as indicated in each 
panel) is plotted versus ip for kF= 200 (solid), 250 (dashed) and 300 (dot-dashed) MeV/c. 
The momentum transfer is q—1 GeV/c. 
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Figure 9: Transverse MEG contribution (solid) compared with the non-relativistic calculation 
of [15] (dashed). 
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propagators: it appears that the dynamical propagator affects the pion-in-ffight more 
than the seagull term (it increases the latter by more than a factor 2 at q=2 GeV/c); 
however, the two effects tend to cancel, so that their net effect is not very significant. 

5.2 Vertex correlations 

Where the v.c. are concerned we have found the following: 

• Their action, while substantial in both the longitudinal and transverse channel, is 
actually dominant in the former by roughly a factor of 3:1, as is evident from Figs. 11 
and 12, where the v.c. contribution (dashed) displayed for several values of q together 
with the free responses (solid). This outcome relates to the minor role played by the 
isoscalar contribution in the transverse response, in turn due to the smallness of the 
isoscalar magnetic moment. Worth pointing out is the oscillatory behavior versus uj of 
the vertex correlations, which induces a hardening of the responses. 

In addition the seagull and vertex correlations tend to cancel in the transverse channel, 
especially for low values of g, whereas for higher q the MEG dominate (see Fig. 13). 
Note that both the seagull and v.c. exactly vanish at the same value of uj, the latter 
coinciding with the QEP for high momentum transfers. 

• In Fig. 14 we study the evolution with q of the v.c. in the longitudinal and transverse 
channel respectively plotting the ratios p^^^^ = R^^p/R^l^} at g=0.5, 1, 1.5, 2, 2.5 
and 3 GeV/c. Glearly the v.c. do not saturate quite as rapidly as the MEG, although 
their behaviour is rather similar and saturation again occurs somewhere above g = 1- 
1.5 GeV/c: thus again scaling of the first kind is achieved at high momentum 
transfers for these contributions. Moreover, similarly to the MEG case, for high q 
the v.c. almost vanish around the QEP (■^^O). In fact the v.c. contributions and the 
seagull MEG contributions are rather similar in shape, but opposite in sign (cf. Figs. 7a 
and 14 (right panel)). Thus, they tend to cancel and, since the relative magnitude of 
the seagull to the v.c. is roughly 2:1, the net effect of the two together is to behave 
the way the former does but with roughly 1/2 its strength. 

• In Fig. 15 we explore the kp dependence of the vertex correlations. They are found 
to grow with kp, much as the MEG do. Again from a semi-relativistic point of view, 
expanding only in powers of hole three-momenta over nucleon mass, we find a leading 
behaviour for kp small that goes as k'jp. In Fig. 15 the ratios and are plotted 
versus ip at q=l GeV/c for kp varying from 200 to 300 MeV/c and we clearly see that 
the v.c. contribution grows with kp. As in the MEG seagull case, for Fermi momenta 
in the typical range (200-250 MeV/c) a slightly lower power of kp'^ actually provides 
the best fit to the full results. The range seen in Fig. 8c reflects the (weak) dependence 
on kp of terms ignored in the semi-relativistic approximation. The basic conclusion 
is similar to that made above for the seagull contribution and hence the total MEG 
at high g, namely, scaling of the second kind is badly broken by effects that go 
roughly as kp. 
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a;(MeV) a;(MeV) 

Figure 10: MEC contribution to versus u with dynamic (solid curves) and static (dashed 
curves) pion propagator at q= 0.5 (a) and 2 (b) GeV/c. The separate pion-in-flight and 
seagull contributions are displayed. 




Figure 11: Correlation current contributions to versus u at g=0.5 (a), 1 (b), 2 (c) and 3 
(d) GeV/c. Solid: free; dot-dashed: self-energy; dashed: vertex correlations. 
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Figure 12: Correlation current contributions to versus u; at q—0.5 (a), 1 (b), 2 (c) and 3 
(d) GeV/c. Solid: free; dot-dashed: self-energy; dashed: vertex correlations. 
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Figure 13: Seagull (solid) and v.c. (dashed) contributions to plotted versus u at q—0.5 
(a), 1 (b), 2 (c) and 3 (d) GeV/c. 
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Figure 14: The ratio (v.c.)/free in the longitudinal (left panel) and transverse (right panel) 
channel plotted versus t/j at kp—^S? MeV/c for various q. The curves are labelled as in 
Fig. 7a. 




Figure 15: The ratio (v.c.)/free in the longitudinal (left panel) and transverse (right panel) 
channel plotted versus ip at q=l GeV/c for various kp. The curves are labelled as in Fig. 7b. 
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5.3 Self energy 



We have found that the self-energy contribution results from a quite delicate cancellation 
between the responses having only the particle or only the hole dressed, as shown in Fig. 16 
in the longitudinal channel for g=0.5, 1, 2 and 3 GeV/c (similar results hold in the trans- 
verse case). This was already pointed out in within the framework of a treatment in 



which relativistic effects were partially incorporated, and it is now confirmed here in a fully 
relativistic context. 

Whereas this cancellation is very substantial at g=0.5 GeV/c, as the momentum transfer 
increases the imbalance between the two contributions grows. Indeed the response associated 
with the particle self-energy is suppressed by the form factors and by the pion propagator, 
but that coming from the hole self-energy is not. As a result, for q > 2 GeV/c the total self- 
energy response almost coincides with the hole result alone and induces a moderate softening 
to the free response. Note that the s.e. contribution does not vanish on the borders of the 
response region. Moreover for high values of u (close to the upper border) it becomes very 
large (see Figs. 11 and 12) and yields a vanishing (or even negative) response, clearly pointing 
to the insufficiency of a first-order perturbative treatment in this kinematical region. This 
effect was already present in the partially relativized analysis of [20| and, significantly, it is 
emphasized by our fully relativistic calculation. Therefore the summation of the full Fock 
series becomes necessary in this case. 

Finally we explore the scaling and superscaling properties of the self-energy correlations 
by studying their momentum transfer- and density-dependence. In Fig. 17 we plot the ratio 

= Rg ^ jRfj.^^ versus ip for the particle (panel a), hole (panel b) and total (panel c) self- 
energies at g=0.5, 1, 2 and 3 GeV/c. As expected, the particle contribution decreases with q, 
going to zero at g ^ 2 GeV/c, whereas the hole contribution, although also decreasing with 
q when not too high, saturates for q >1 GeV/c. As a result the total self-energy, displayed 
in panel (c), grows with q (in contrast to all other cases considered above) in the range 
g=0.5-2 GeV/c, then stabilizes typically at about 30-40% of the free response to the left of 
the QEP, thus inducing the above-mentioned softening of the longitudinal response. Similar 
results are found in the transverse channel. In summary, again scaling of the first kind 
is achieved at momentum transfers somewhat below 2 GeV/c. 

In Fig. 18 we display the same ratio for three different values of the Fermi momentum at 
q=l GeV/c. It appears that, as for the other two-body correlations, the self-energy relative 
contribution grows with kp, specifically i?fg oc kp, although not uniformly in ip (see Fig. 8, 
panel d) — recall that in the first-order analysis presented in this paper the edges of the 
response region are not treated adequately for the self-energy contribution and thus should 
not be taken too seriously. Where the self-energy contribution is correctly modeled (away 
from the edges) we again see breaking of second-kind scaling by roughly kp. 



To complete the presentation of our results we display in Figs. 19 and 20 the global 
responses in first order of perturbation theory and compare them with the zeroth order 
ones (free responses) for several momentum transfers. Here one assesses the impact of the 
global two-body current contribution to the responses. First the overall effect of the two- 
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Figure 16: Particle (dashed) and hole (dot-dashed) contributions to the longitudinal 
self-energy (solid) versus lu ior q = 0.5 (a), 1 (b), 1.5 (c) and 2 (d) GeV/c. 
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Figure 17: The ratio of the particle (a), hole (b) and total (c) s.e. contribution to the free 
longitudinal response is plotted versus for q—0.5 (solid), 1 (dashed), 2 (dot-dashed) and 3 
(short-dashed) GeV/c. 
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Figure 18: The ratio of the particle (a), hole (b) and total (c) s.e. contribution to the 
free longitudinal response is plotted versus -0 for A;i;'=200 (solid), 250 (dashed) and 300 
(dot-dashed) MeV/c and q=l GeV/c. 
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Figure 19: Longitudinal response versus cu including all first-order contributions (dashed) 
compared with the free result (sohd) at q—0.5 (a), 1 (b), 2 (c) and 3 (d) GeV/c. 
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Figure 20: Transverse response versus uj including all first-order contributions (dashed) 
compared with the free result (sohd) at g=0.5 (a), 1 (b), 2 (c) and 3 (d) GeV/c. 

body currents appears sufficiently modest to justify our first-order treatment. Next the 
softening at large q appears to be common to both L and T channels, whereas at low q the 
longitudinal response displays a hardening that is absent in the transverse one. Also evident 
is the already-noted almost vanishing of the two-body correlation contribution at the peak 
of the free responses. Finally the unrealistic dominance of the self-energy contribution on 
the upper border is apparent. 

6 Conclusions 

In this paper we have addressed some of the ingredients that enter when attempting to 
model the longitudinal and transverse nuclear response functions for inclusive quasielastic 
electron scattering. Because of its complexity this problem cannot presently be viewed as 
solved: indeed, although many papers claim success in accounting for the data, a closer 
scrutiny reveals a different situation. It is not only that contributions left out in various 
analyses are far from being small, but, even more serious, fundamental physics principles 
(Lorentz covariance, gauge invariance and unitarity) turn out patently to be violated. Thus 
the successes in reproducing the experiments, regrettably few both in connection with the 
L/T separation and with the set of nuclei or range of q explored, often reflect more an 
adjusting of parameters than a real understanding of the physics involved in the QEP. 

Searching for an improvement of this situation we have engaged in a more systematic 
approach to the problem, assuming as a zeroth-order approximation the RFG model. This 
appears justified on two counts: i) because RFG is Lorentz covariant, ii) because the surface 
of the nucleus is largely irrelevant in the QEP. 

One of our basic themes has been the implementation of Lorentz covariance. In past 
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work |[T^ we attempted to approximate the full theory by identifying a dimensionless variable 
that is small enough to be suitable to use in setting up a non-relativistic expansion of the 
responses (namely the momentum of a nucleon lying below the Fermi surface compared with 
its mass). In contrast, in the present work no expansion whatsoever is involved and our 
calculation fully respects relativity. The motivation for this should be clear: the theoretical 
analyses of the present and future experiments at TJNAF demand as full an implementation 
of relativity as is feasible. 

Since the physics of the QEP is remote not only from the physical surface of the nucleus, 
but from the Fermi surface as well, we consider a treatment in terms of nucleonic and mesonic 
degrees of freedom (the latter viewed both as force and current carriers) to be appropriate. 
Hence, as a first step, we focus on pions as they can be expected under these conditions 
to be an important, if not the major, carrier of the currents that respond to an external 
electromagnetic field impinging on the nucleus. The pions, in our framework, are dealt with 
in first order of perturbation theory, since their effects on the free responses of the RFG are 
not expected to be too disruptive. 

A point worth making in connection with all of the two-body currents considered in the 
present work is the following: their contributions appear to be quite small (if not vanishing) 
at the peak of the quasielastic responses. This corresponds to viewing the Fermi sphere to 
be split into two domains, an inner one associated with a sphere with a radius smaller than 
kp, and the other corresponding to the spherical shell representing the difference between 
the two concentric spheres. The contributions tend to cancel one another and at the peak, 
where the whole Fermi sphere is involved in the response, they nearly do so completely. Of 
course the radius of the inner sphere and the sign of the contributions are current dependent. 

Gauge invariance is a fundamental property we have also addressed in this work. In 
particular we have verified whether or not the continuity equation it implies is valid order 
by order in perturbation theory. Notably we have succeeded in showing that the continuity 
equation for the one-body (single-nucleon) and the two-body (MEG and correlations) cur- 
rents is indeed fulfilled; hence our approach deals consistently with forces and currents. At 
the same time we have also shown how crucial relativity is in achieving this and how much 
care is needed in dealing with the ph matrix element of the self-energy current. Actually the 
present approach points to the necessity of renormalizing both the energies and the spinors 
in the self-energy ph matrix elements. Work under these lines is under progress p3 |. 

As far as the self-energy contribution is concerned, the present results show that for high 
momentum transfers a first-order perturbative treatment seems to be insufficient to deal 
with the kinematical region corresponding to high values of u (close to the upper border of 
the response region). This shortcoming demands the inclusion of higher-order perturbative 
terms, namely the summation of the whole Fock series, a task performed in a partially 
relativized framework pO| but not yet achieved in a fully relativistic one. It might be worth 
pointing out that the above shortcoming relates to the "wrong" analytic properties of our 
first-order polarization propagator: it should be a meromorphic function of the energy, yet 
it displays a double pole. 

With regard to the other contributions of the two-body currents to and R^, our 
results are highly satisfactory. Indeed for the MEG we have found that their contribution 
is small enough to be well handled in first order. In particular both the pion-in-fiight and 
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seagull contributions are very small in the L channel where the virtual photon exchanged 
between the electron and the Fermi gas couples to the charge of the pion. We thus see 
that the MEC affect the Coulomb sum rule very marginally. In the T channel instead the 
MEC are more significant and the following features emerge (see also Section |^ for itemized 
observations): 

• the seagull contribution dominates, in accord with the spin nature of the MEC physics, 
and displays an oscillating behaviour versus ui yielding a (modest) softening of the 
response; 

• the MEC contribution does not vanish when q increases, but ultimately at very high q 
it displays scaling of the first kind; 

• the MEC contributions break the scaling of the second kind; i.e., they go as higher 
powers of kp than does the RFC. In fact, the total MEC contribution goes roughly as 
k\ with respect to the RFC and thus when kp becomes small the MEC contribution 
tends to vanish. Hence, in going from the lightest nuclei where kp is as small as 
about 55 MeV/c (viz., for the deuteron) to heavy nuclei where it may be as large as 
260 MeV/c, one should see the Ip-lh quasielastic MEC contribution go from being 
essentially negligible to becoming typically 5-10% of the total transverse response. 

With regard to the correlation contribution arising from the vertex corrections (v.c), we 
have found that the L channel dominates over the T channel by a margin of roughly 3:1. 
The longitudinal response effectively picks up only these correlation contributions, since the 
MEC effects are so small there, and the former contribute to the total at roughly the 10- 
15% level. Indeed, were these to be the only contributions needed in addition to the RFC 
response itself, then we would expect the total to shift in uj (or in ip). Since the correlation 
contributions are roughly symmetrical about the QEP, it is important to note that their 
impact on the Coulomb sum rule should be very small, perhaps only at the few percent 
level. 

The correlation contribution to i?^ is similar to the MEC contribution, but is smaller, 
roughly 1/2 the size of the latter. Since the two are of opposite sign they tend to cancel 
and thus the total is similar to the MEC contribution but is cut down by a factor of two. 
This statement is of course only approximate, since the correlation contribution attains its 
asymptotic value as a function of q somewhat less rapidly than does the MEC contribution. 
Interestingly, both the seagull and v.c. contributions appear to vanish at exactly the same 
value of UJ. Overall the total (the sum of Ip-lh MEC + Ip-lh correlations) yields contribu- 
tions to be added to the RFC that (1) do not go away as q becomes very large, (2) attain 
scaling of the first kind at momentum transfers somewhat above 1 GeV/c, while (3) they 
clearly violate scaling of the second kind, going approximately as k\ relative to the RFC 
result (which of course scales). 

In conclusion, in this paper we have pursued the program of a systematic analysis of the 
inclusive nuclear responses in the domain of the QEP and beyond. The motivation for this 
effort lies in the recognition that in this energy domain, unlike for the ground and low-lying 
excited states of nuclei, the interaction among the constituents of the system cannot be 
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subsumed into a potential. Currents are important; and moreover the consistent treatment 
of currents and forces is important. For the first time one is able to investigate the scaling 
behavior of at least some of the contributions to the inclusive cross section under high-energy 
conditions using a non-trivial (interacting) model in which issues relating to covariance and 
gauge invariance are addressed consistently. 
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Appendix A. The gauge invariance of the two-body currents ph 
matrix elements 

Following the study presented in Section |^ of gauge invariance at the level of the free- 
space particle-particle matrix elements, here we extend the analysis to the particle-hole 
channel, deriving the contribution to the continuity equation of the isoscalar and isovector 
s.e., v.c. and MEC particle-hole matrix elements. We start by evaluating the divergence of 
the correlation particle-hole matrix element {ph~^\j'^\F) for the s.e. and v.c. contributions; 
next we address the MEC ph matrix elements. 

• Self energy 



From eqs. (25,26[) we get 



O 7i - ^ (f-^)(^-m)(/>-^) 

2mV^mi JEpEh k [P - Ky - 



O V - m m (^- fl){^-m){^- ft) 



Note that Fi cannot be taken out of the matrix element since it acts on the isospinors. Now 
from the relations 

Sf{P) Qu{h) = u{h) (58) 

w(p) ^Sp{H) = -n(p) (59) 

m(p)(/'- ^)(^-m) = 2mu{p){f- ^) (60) 

-m){^- /I)u{h.) = -2m{^ -m)u{h) (61) 

the following expressions are derived: 

^ 3P m ^ m , AK-m)iP-K)^ , 

V^mi JEpEh k [P-Ky -mi 
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Q-nh = --—^^== 2^ — M p Fi— — — h . 63 



Vertex correlations 



From eqs. ( P7| , pBD the four-divergence of the v.c. matrix element is found to be 

V^mi JEpEh li^ {K - Hy - mi 



V^mi JEpEi, k [P-Ky -mi 



We now exploit the identities 



Sp{K + Q) Q{^ + m) = +(^ + m) (64) 
+ m) iQSpiK -Q) = -i^ + m) (65) 



to get finally 



n T rn m (^- P){^-m) 

r^m^ y'EpEh k [K - Hy - mi 



^ ^ f m ^ m ^ ^ ^ (K-m)(P- ^) „ , 
^^"^^ /EpEh k ^k (^ - /'^)^ - 



If the expressions (^,^,^,0) are split into their isoscalar and isovector parts, as illus- 
trated in Section we get 

n W(^) - (5) m m (^-myp- ^) 

"^'^^ - -y^f' -T^E^-lP) ^p^Ky-mi 



n -wCV) _ 3/^ m (^-m)(/^-^)r3 

P F^m^ ^EpEi, k ^k (^^ - ^)^ - mi 

'^■'''^ - "^"^^ 71^?^"^^^ (ir-i/)^-< 

V"^m^ y'EpEh k ^k (^ - Hy - mi 

l^^m^ JEpEh k ^k (^ - ^^)^ - mi 



34 



Q ■ 

Q ■ fi(^) 



3f^ /Q\ m 



P 



2^ T^«(p)^75 — —9. — z:^«(h) 



iV) 



m 



EpE^ r 



(P - - ml 
(^-m)(/>- ^) 

{p-xy-mi 



From these relations we learn that: 

• In the isoscalar channel the self-energy and vertex contributions cancel 

Q ■ ) + Q ■ = Q . T^f ) + Q . ^(^) = 0. 



Ts + i£3aferarb)M(h). 



(6^ 



This is in contrast with the non-relativistic result where the self-energy is by itself 
gauge invariant. 



In the isovector channel we get 

2p 



2P 



2^ T^"(P) 775 — «W 



^P^h k ^k 



'p 
m 



vat, 



ah 



^p^h r ^k 



These expressions, using the Dirac equations /fu(h) = mM(h) and m(p) /' = mu(p), can be 
further simplified to yield the following four-divergence of the correlation current 

1 



2f 



P 

m 



E 



EpEh k ^k 



K -P-m^ 



K-H -m'^ 



{P-Kf-ml {K-Hf 



mt 



nuihi). 



This contribution is exactly cancelled by that of the MEC (seagull and pion-in-flight) as we 
illustrate in what follows. 

• MEC 

Using the expressions given in Section |3.2| for the ph matrix elements corresponding to the 
seagull and pion-in-flight currents in eqs. ( pl] , |22D , the associated four-divergences are found 
to be 



-^1 ^£3, 



m 



•ah 



X 



-Sp-^h 



+ 



— m) 



k 

2mp 

m 



Kf -ml {K- Hf - ml 



M(h) 



EpEh 

{Q^ + 2H ■Q-2K ■ Q) 



^ E^\{P-Kf-ml][{K-HY-ml] 



u{p)Ta{^ - m)TbU(h). 
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Exploiting the Dirac equation and after some algebra the above can be recast as follows 

{Q ■ ?^^^) = {Q ■ r) + {Q ■ F) 

2/2 m „ m , f K -P-m^ K ■ H - 1 

Fi 'tesab I = 2^ —u[T£>)Ta \ — — ^ ^ - — — ^ ^ \ nu[h.) . 



We have thus proven that the correlation and MEC ph matrix elements satisfy current 
conservation, i.e., {Q ■ j^) + {Q ■ j^'^^'^) = 0. 

Appendix B. Poles in the vertex correlation diagrams 

• Forward Diagram 

As shown in Section p.2|, the forward-going diagram (Fig. 2d) of the v.c. ph matrix element 
involves integrals of the type 

The question then is if, for given q and u, a four-momentum K^^, with k < kp, exists such 
that P' = K + Q is on-shell. The three-momentum k lies on the sphere A with radius k (see 
Fig. 21). On the other hand the three-momentum p' should lie on the sphere C of radius 



p' = v/(^k + ^)2-m2 , (70) 

where E'k = y/k"^ + mP' is the energy of the hole. Moreover, since p' = k + q, it should also 
lie on the sphere B which corresponds to A translated by an amount q. The pole is clearly 
hit when the two spheres B and C intersect. 

A quantitative assessment of the position of the poles can be more easily derived in the 
non-relativistic framework. In this case indeed the particle momentum reads 

p'2 = fc2 ^ 2muj . (71) 

Let us now keep q and k fixed and study the behaviour of the poles as a function of uj. On 
the basis of eq. ( |7TD we can consider three situations: 

• i) UJ = q^ /2m (QEP). Here p'^ = k"^ + q^ and consequently k is perpendicular to q, so 
the pole appears for cos6'fc = 0. This situation is illustrated in Fig. 21(b). 

• ii) > q^ /2m (above QEP). The particle momentum satisfies p'^ > k'^ + q^ (Fig. 21(a)). 
In this situation the radius of the sphere C is larger than in the previous case (Fig. 21(a)), 
and therefore 6^ < ^ and cos6'fc > 0. 

• iii) UJ < q^/2m (below QEP), i.e., p'^ < k^+q^ (Fig. 21(c)). Here Ok > ^ and cos^^ < 0. 
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(a) 



(b) 





Figure 21: Poles in vertex correlation. Diagrams (a-c) refer to the forward-going situation 
and (d) to the backward-going one. 



To estimate the minimum value of the hole momentum k allowing for the existence of a 
pole, from Figs. 21(a) and 21(c) it clearly appears that, for given q and u, there are no poles 
if the value of k is very small: in fact in this case the spheres B and C do not intersect. 
More precisely, the range of allowed values of k is limited by the condition (for k < q): 



q — k<p<q + k. 
Within the non-relativistic framework (see eq. ([7T|) ), this condition reads 

|mAci;| < qk , 



(72) 
(73) 



where Auj is defined as the cu-shift from the QEP: Auj = on — Then the minumum value 
of k allowing for the existence of a pole is 



kr, 



m 

q 



q' 




A 


UJ 


= m 


K 


2m 




K 



(74) 
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It then follows that at the QEP one has fcmin = and therefore the pole always comes into 
play. 

• Backward Diagram 
In the case of the backward-going v.c. (Fig. 21(d)) the integral to be considered is 

J [J^ — Qj — m'^ + le 

Also here one can study the conditions that allow a particle with four-momenutm P' = K—Q 
to be on-shell for given q and uj. Now clearly the energy associated with p' is E'^ = — uj 
with E\^ = y/k'^ + m^. The momentum p' must lie on the sphere C of radius 

p' = ^J{Ei, - cj)2 - m2 <k<kF , (76) 

which now is contained inside the sphere A (see Fig. 21(d)). On the other hand, since 
p' = k — q, it must also lie on the sphere B obtained for A by translating by an amount — q. 
For q > 2kp (no Pauli-blo eking) the intersection between B and C vanishes: hence no poles 
can be found. 



Appendix C. Polarization propagator with nucleon self-energy 

In this appendix we derive the general expression for the self-energy polarization propaga- 
tor with any number of insertions, S(P), in the particle and/or in the hole lines. From 
the obtained expression one can derive, as particular cases, the leading-order response (no 
interaction lines) and the first-order self-energy response (with one interaction line). 

Let us then consider the polarization propagator containing n self-energy insertions S(i7) 
in the hole line and / insertions S(P) in the particle line. It is given by 



n 



-iTi 



dhf^d^h 



T^'{Q)[S{H)T.{H)YS{H)T'{-Q)[S{H+Q)T.{H+Q)]'S{H+Q). (77) 



Using the nucleon propagator in the medium in eq. (^2|) we can express the product of n -|- 1 
propagators appearing in eq. (|77D as a derivative of order n according to 



[S{H)T.{H)fS{H) = 

= [(/r + m)S(i7)]"(/r + m) 
= [{/I + m)^{H)]"iP + m) 



e{h - kp) 



+ 



9{kF - h) 



1 



a=0 



e{h - kp) 



+ 



e{kF - h) 



a 



m? + ie H"^ — a — rn? — ieho 



1 d" 



n\ doi"- 



a=0 



a 



+ ie 



+ 2ni9{kF - h)6{H^ -a- m^)d{ho) 



(7J 
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A similar equation holds for the propagation of a particle. The global polarization propaga- 
tor, with both the particle and the hole line dressed, can then be written as 



n 



1 (i" 



nW. da 
X Tr 



a=0 

dhod^h 



dpi 



(3=0 



X 



(27r)4 
1 



r^(g)[(^ + m)i:{H)f{p + m)r'^(-g)[(/' + m)T.{P)]\p + m) 



— a — m? + ie 
1 

P"^ — [3 — rn? + ie 



+ 2-Kie{kF - h)5{H^ -a- m'^)9{ho 
+ 27rie{kF - p)S{P'^ - ^ - m'^)e{po) 



(79) 



where P = H + Q. 

It gives rise to four terms. One of these contains the product of the two free propagators, 
namely 

I \ f I \ 

(80) 



— a — + it) \P'^ — (3 — rn^ + ie J ' 

and yields the divergent vacuum contribution, which should be subracted out and pertains 
to a domain beyond nuclear physics. Then, after some algebra, one obtains 



27r 









da"" 




I 

13=0 •' 



dhod^h 

'WW 



aH,P,Q) 



e{kF - p)6{P^ -13- m^)e{pol ^ e{kF - h)6{H^ -a- m^)e{ho) 



if ^ — a — + ie if ^ — (3 — ''^^ + 

+ 2nie{kF - p)e{kF - h)5{P'^ - f3 - m^)5{H'^ -a- 711^)6 (Pq) 6 {Hq 



where the function 



^niiH, P, Q) = li^iiho, h, po, p; go, q 
1 



nW. 



rTr 



r^(g)[(^ + m)s(ii)]"(/r + m)r"(-g)[(/' + m)s(p)]'(/' + m) 



(82) 



has been introduced. According to eq. (H) we then obtain for the nuclear tensor 



— im [n'^'^ - n(; 

TT 



27r 





d' 




da" 


a=o d(3^ 


I 

13=0 



dhod^h 



liiiiH, P, Q)6iP' -13- m')5{H' - a - m') 



X [e{kF - p)e{po) + 9{kF - h)e{ho) - 2e{kF - p)e{kF - h)e{po)e{ho)] 
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Exploiting now the identity 



[e{kF - p)e{po) + 9{kF - h)9{ho) - 29{kF - p)e{kF - h)9{po)9{ho)] 
9{kF - h)9{ho) [1 -9{kF- p)9{po)] + 9{kF- p)9{po) [1 - 9{kF - h)9{ho)] 



we can recast the above according to 



1 



- -Im [n'^^ - Ho 



^5) 



where 
and 



27r 





d' 




da"" 




I 

/3=0 



dhnd^h 



CiH, P, Q)6{P^ -(3- m^)6{H^ - a - m^) 



(2vr)4 

X 9{kF - h)9{ho)[l - 9{kF - p)9{po)] 



(86) 



2n 









da"" 




I 

(3=0 



dhnd^h 



(if, P, Q)5(P' -P- m^)6iH^ -a-m" 



(2vr)^ 

X 9ikF - p)9ipo)[l - 9ikF - h)9iho)] 



^7) 



Now eq. ( p6D corresponds to the usual hadronic tensor for electron scattering, whereas eq. (p7| ) 
corresponds to a process with negative energy and momentum transfer, hence it should be 
disregarded. 

Next the integration with respect to can be explicitly performed by using the second 
of the 5-functions appearing in eq. (|5B|). After some algebra one then gets for the hadronic 
tensor the expression 





d^ 




f d'h /^r(^h(«),h;S(,(a) + go,p;g) 




„=o dP^ 


13=0^ 





2E' {(3) " ~ ' ^ ^ 

where 

E{^{a) = Vh'^ + a + m? (89) 

E'^{(3) = ^p' + P + m?. (90) 

Finally, exploiting the first (5-function in eq. (|86D , one can cast the response function in 
the form 





d' 




f d'h /„T(^h(«) 


,h;E'^{P),p;q) 




n dp^ 


13=0 


(27r)3 AE', 


»Ei,{P) 



x6iE^{a) + go - E'^iP^kF - h)9{p - kF). (91) 
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For n = I = one gets the well known free hadronic tensor 



= / + ^0 - E^)eik, - h)eip - k,) , (92) 



for n = 1, Z = the first-order hole self-energy hadronic tensor 



da 



/ JTT^ ^n.// ^i^hW + ?o - Ep)e{kF - h)e[p - kp) (93) 

a=o-^ (27r)-^ AE^[a)Ep 



and finally for n = 0, / = 1 the first-order particle self-energy hadronic tensor 
d^h JoT(i?h,h;E;,(/3),p;g) 



dp 



(2^ AE^) 5iE^ + %-Epm0ik,-h)eip-k,). (94) 



Note that the last two expressions coincide, as they should, with the ones already given in 
Section |[ 



The hadronic tensor in eq. (|91|) is one of the {n + /)-th order contributions to the full 
Hartree-Fock polarization propagator, which should be obtained as an infinite sum of all 
perturbative orders. It should be stressed that in a basically non-relativistic framework the 
HF series can be summed (see, e.g., [^), whereas in our fully relativistic framework the 
summation is far more difficult to perform because of the matrix nature of the self-energy 



S |6|. 
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